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ABSTRACT 


A  theory  is  presented  for  an  end  effect  in  the  current  response 
of  a  highly  negative,  cylindrical  Langmuir  probe  in  a  collisionless 
plasma  flow.  Under  conditions  where  the  ratio  of  probe  radius  to 
Debye  length  is  small  and  the  ion-acoustic  Mach  number  is  large,  the 
current  exhibits  a  strong  peak  when  the  probe  axis  is  brought  into 
alignment  with  the  flow  direction.  Closed  formulae  are  given  for 
the  maximum  and  angular  half-width  of  the  peak,  and  universal 
graphical  results  are  presented  for  the  entire  peak  structure.  The 
theory  shows  very  good  agreement  with  experimental  data.  The  use  of 
the  end  effect  for  diagnostic  purpose?,  in  particular  for  the  deter¬ 
mination  of  the  ion  temperature,  '.s  discussed. 
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I  -  Introduction 

The  theory  of  an  Infinitely  long,  cylindrical  Langmuir  probe  in 

a  collisionless,  quiescent  plasma  is  particularly  simple  when  the 

probe  radius  r^  is  smaller  than  the  Debye  length  X^,  since  then  the 

1  2 

current  is  "orbital  motion  limited"  and  the  old  Langmuir  analysis 

is  valid.  Actually,  that  analysis  is  valid  even  in  the  more  general 

case  in  which  the  plasma  is  in  motion  relative  to  the  probe.  If  0 

is  the  angle  between  probe  axis  and  flow  direction,  and  -eV  and 

P 

2 

nijU  »  kT^,  (where  is  the  applied  potential,  icTe  is  the 

2 

electron  thermal  energy  and  kT^  and  m^U  are  the  ion  thermal  and 
directed  energies,  respectively)  the  current,  J  ,  to  a  probe  of 
length  i,  as  l  -*■  <*»,  is  given  by 


J  ll  -  2N_eU  sin  0  r  (1  -  2Z.  eV  /mV  sin2  0]1/2  (1) 

“  U  p  i  p  i 


vhere  is  the  plasma  density  and  Z^  is  the  ion  charge  number. 

Notice  that  according  to  Eq.  (1)  decreases  monotonically  as  0 
goes  from  tt/2  to  zero. 

In  any  actual  experiment,  however,  1  must  be  finite.  Recently, 

3 

probe  current  data  have  been  reported  from  both  satellite  and 
4  2 

laboratory  experiments,  for  rp  «  XD,  m^U  »  and  Jl/rp  as  high 

as  820,  that  show  a  striking  disagreement  with  as  given  in  Eq.  (1). 
The  current  observed,  J,  was  close  to  Jw  as  long  as  0  was  not  small, 
but  as  the  probe  approached  the  aligned  orientation  J  exhibited  a  sharp 
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rise  that  peaked  at  8  *  0  at  a  value  many  tines  larger  than  J##(6  •  0) . 

3  5 

This  phenomenon  may  be  explained  *  as  an  end  effect  due  to  the  finite 
length  of  the  probe.  This  shows  that  extremely  long  probes  may  be 
necessary  if  Eq.  (1)  is  to  be  applicable  in  the  interpretation  of 
probe  characteristics. 

Of  more  interest  for  diagnostic  purposes,  however,  is  the  end 
effect  in  itself.  The  peak  may  be  quite  strong,  and  should  be  possible 
to  use  it  in  the  determination  of  the  relative  direction  of  the  plasma 
flow,  and  of  a  number  of  plasma  parameters.  Of  particular  interest 
is  the  fact  that  both  the  height  and  the  half-width  of  the  peak  are 
often  sensitive  to  the  ion  temperature.  This  is  very  Important  be¬ 
cause  no  other  feature  of  probe  response  is  known  to  be  noticeably 
dependent  on  T^. 

Bettinger  and  Chen  were  the  first  authors  to  present  a  theoretical, 

although  crude,  analysis  of  the  end  effect;  an  important  limitation 

of  their  approach,  as  pointed  out  in  Ref.  5,  was  that  l  had  to  exceed 

2  1/2 

a  minimum  value  l  [1  a»  3A_  (m.U  /Z.icT  )  ].  For  i.  <  l  some 

m'm^D  i  ie  m 

numerical  computations  were  carried  out  by  Hester  and  Sonin  for  0  »  0. 

The  present  analysis  starts  from  a  similarity,  suggested  in  Ref. 

5,  between  the  present  steady-flow  problem  and  a  time-dependent  one 
involving  a  quiescent  plasma.  In  the  next  section  the  conditions  for 
the  validity  of  this  Hester-Sonin  (H-S)  similarity  are  discussed  in 
detail  and  the  basic  points  of  a  theory**  recently  developed  for  the 
time-dependent  problem  are  Introduced;  this  theory  illustrates  clearly 
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Che  anomalous  behavior  of  the  current  as  a  function  of  8.  In  Sec.  Ill 
analytical  and  graphical  results  for  the  main  features  of  the  end  effect 
are  presented  and  compared  with  experimental  data.  The  applications 
of  the  effect  are  discussed  In  Sec.  IV,  and  Bettlnger  and  Chen's 
analysis  is  discussed  In  an  Appendix. 

II  -  Basic  Formulation 


We  consider  a  long,  cylindrical  Langmuir  probe  with  length  l  and 
radius  rp  In  a  collisionless  plasma  with  unperturbed  thermal  energies 
icTe  and  kT^,  density  and  bulk  velocity  relative  to  the  probe  U  at 
an  angle  8  with  Its  axis.  The  probe  potential  Vp  Is  negative  and  such 
that 


m  U2  <<  icT  <<  -  eV  ;  (2) 

e  e  p 

the  electron  current  is  then  negligible  and  the  perturbed  electron 
density  Is  given  by  Boltzmann's  law 

Ne  »  Nq  exp  (-  <j/)  (3) 

where  <J>  =  -  eV/icTe  is  the  nondimensional  potential  field.  Defining 

0  -  We  >  M  "  U/Ad Wp  ,  (4) 

l  -  t/XD  ,  e  -  rp/XD  ,  (5) 


■'  -  ■‘-‘-r-.v  r 
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we  assume  that  M  and  l  are  large,  c  small  and  6  S  1;  M  is  the  lon-acoustic 

2  1/2 

Mach  number,  A^  =  (tc^/AirN^e  )  is  the  electron  Debye  length, 

2  1/2 

Up  =  (4itNqZ^6  /m^)  is  the  ion  plasma  frequency  and  m^  and  Z^  are 
the  ion  mass  and  charge  number.  If  j(z)  is  the  current  density  at 
the  probe  at  distance  z  from  its  tip,  the  average  current  density 

i 

j  ■  i"1  j  j(z)  dz 
0 

can  then  be  written  as  a  nondimens ional  function 


4“  “  4“  U,  M,  e,  0,  *  ,  0) 

J.  P 


(6) 


where  j  =  j  (1  -*■  ») ;  according  to  Eq.  (1) 


j,,,  -  ir"1  NQe  U  sin  0[1  -  2Z±  eV  /z^U2  sin2  0]1/2, 


(7) 


The  total  current  to  the  probe  is  J  ■  2irr  i  jM  times  j/j^.  To  obtain 
this  last  quantity,  the  ion  Vlasov  equation  must  be  solved  together  with 
Eq.  (3)  and  Poisson's  equation.  The  last  one  reads 


1  3_ 
P  9p 


P 


14  +  I_  AJfi  +  Li  iJfc 

3P  p23*2  i23?2 


2 

e  v 


(8) 


we  have  introduced  cylindrical  coordinates  r,  4>  and  z  and  have  defined 
dimensionless  quantities 


-5- 

p  -  r/rp  ,  C  -  z/l  ,  \>±  »  (Z±N1  -  Ne)/NQ  (9) 

where  is  the  perturbed  ion  density.  equals  ^  =  -  eVp/xTe  at 

the  probe  and  zero  at  infinity.  The  ion  distribution  function  far 

1/2 

ahead  of  the  probe  must  also  be  known  [even  though  U  >>  (xT^/m^)  ] 

and  here  will  be  assumed  to  be  Maxwellian.  [Thermal  velocities  may, 
and  will,  be  neglected  in  the  motion  along  the  z-axis,  but  thermal 
motion  in  the  p  -  <p  plane  is  of  fundamental  importance  when  0  is 
small. ] 

Let  us  begin  by  considering  the  limit  6*0  (which  also  implies 

3/3$  *  0).  Hester  and  Sonln^  studied  this  limit  and  pointed  out 
2  2  2 

that,  if  SL  and  M  /ipp  are  so  large  that  3  i|//3C  can  be  neglected  in 
Eq.  (8)  and  the  ion  velocity  along  the  z-axis  can  be  well  approximated 
by  its  unperturbed  value  U,  the  steady  flow  problem  is  equivalent  to 
-  time-dependent  o^e,  wherein  an  Infinitely  long  probe  is  Immersed  at 
time  t  *  0  in  an  unperturbed  quiescent  plasma,  all  other  conditions 
being  the  same  of  the  original  problem.  The  time  of  flight  of  the  ions 
down  the  probe  z/U  and  the  current  density  at  z  are  respectively  equi¬ 
valent  to  the  time  t  and  the  (spatially  uniform)  current  density  at  t. 
As  z  (t)  Increases,  the  ion  distribution  function  readjusts  Itself  and, 
if  1  is  large  enough,  the  "infinite"  probe  (steady  state)  limiting 
current  density  j  (8  *  0)  will  eventually  be  reached.  In  the  context 
of  the  time-dependent  problem  we  can  write 

i(6-0,--U  (10) 

0 


where 


tt  -  m 


(u) 


is  the  time  equivalent  of  the  length  of  the  probe  in  the  flowing 
plasma. 

The  time-dependent  problem  has  been  recently  analyzed  by 
Sanmartin*’.  His  approach  is  based  on  the  following  points:  (I)  It 
is  possible  to  derive  an  accurate  expression  for  the  electric  field 
3«p/ 3p  (for  limited  values  of  p)  without  solving  simultaneously  the 
ion  Vlasov  eq  tion.  Integrating  Polsscn's  equation 


1  _ 
0  3p 


3^1 

3p 


2 


e  v 


yields  b  formal  expression 


3i|) 

3p 


_ 


6(t) 


+  (P2  -  1)  <v(p,t)> 


(12) 


(13) 


where  -ij; 

P 


6(t)  is  the  field  at  the  probe  at  time  t  and  <v>  is  defined  by 


2 

P 

(p2  -  1)  <v(p,t)>  -  |  d(p'2)  v(p',  t) 

1 


(14) 


5(t)  is  found  to  change  very  little  from  t  ■  0  to  t  *  »  so  that  an 
intermediate  constant  value  6  can  be  used  in  Eq.  (13);  one  finds 


6  1  »  In  e"1  +  Y(c,  i>  ) 

P 


(15) 


mmmmmmHmmtesmHKmmsmm 


where  Y  is  given  in  Fig.  1.  <v>  varies  over  the  entire  range 

1  *  <v>  5  0;  however,  in  a  certain  neighborhood  of  the  probe,  roughly 
—  1/2  -1 

P  $  Pm  =  (2ij»p  6)  e  ,  <v>  is  found  to  change  very  slowly  and  to 

have  a  central  value  v  that  is  approximately  0.80.  For  details  on 

these  results  and  their  accuracy,  see  Ref.  6.  Thus,  for  p  s  p  ,  the 

m 

electric  field  may  be  correctly  approximated  by  a  function  that  only 
depends  on  p 


3jj» 

3p 


i) 


(16) 


Although  for  e  small  the  concept  of  a  sheath  has  little  meaning 
("orbital  motion  limited"  current  implies  an  infinite  sheath  in 

2  — i 

the  Langmuir  sense  )  both  f  and  exp(6  ),  which  are  close  to  each 

other,  may  be  thought  of  as  characteristic  sheath  radii.  If  t  is 

m 

the  typical  time  of  flight  to  the  probe  of  ions  that  were  at  the 

boundary  of  that  sheath  at  t  ■  0,  it  is  clear  that  the  theory  is  only 

valid  for,  roughly,  t  $  t  . 

m 

(II)  With  the  field  known,  ion  trajectories  may  be  computed 
explicitly.  Moreover,  the  current  to  the  probe  is  linear  in  the 
(unperturbed)  ion  distribution  function  at  t  ■  0,  fg(v*).  su^ices 
therefore  to  determine  the  current  for  the  simplest  possible  f^,  that 
for  which  all  ions  have  velocities  of  the  same  magnitude  v*  and 
direction  (parallel  to  s.n  arbitrarily  chosen  polar  axis);  see  Fig.  2. 
Once  that  current,  j*(v*)/j  ,  has  been  found,  the  current  for  any 

A  .  ™ 
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other  fQ  is  given  by  a  definite  integral,  J  fg(v*)  dv*  J *(v*)/ja«  To 
find  J*/ one  can  use  energy  and  angular  momentum  conservation  to 
divide  the  p  -  <J>  plane  at  t  ■  0  in  two  mutually  exclusive  regions 
A*(v*)  and  B*(v*):  a  point  p  -  4>  belongs  to  A*  if  an  ion  having  such 
initial  coordinates,  and  moving  under  the  field  given  by  Eq.  (16), 
will  strike  the  probe  at  some  t  >  0.  The  time  of  flight  to  the  probe 
of  every  point  in  A*  can  be  computed  and,  therefore,  a  function 
a*(v*,  t)  can  be  determined  which  represents  the  area  of  that  part  of 
A*  which  has  been  "collected"  by  the  time  t.  The  current  per  unit 
length  of  probe  is  then 


2irrp  j*(t)  -  NQe  da*/dt 


(17) 


and  defining  a(t)  ■  j  a*  dv*  f0(vf)t  we  have 


2irrp  j(t)  -  NQe  da/dt 


(18) 


The  average  value  of  j  is 


T.V_a(tl 

J  2irr  t 


(19) 


and  in  nondimens ional  form 


i-  -  _ 

(8<j>  )1/2  T 

P 


(20) 
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A  A  2 

where  t  ■  w  t  and  a  Is  a  nondimenslonal  form  of  the  area,  a  *  a/r  . 

P  P 

The  general  behavior  of  j/j^  can  be  determined  by  a  qualitative 
discussion  of  j*/j(D  assuming  that  v*  is  of  the  order  of  the  charac¬ 
teristic  ion  velocity.  The  condition  for  a  point  to  belong 

to  A*  follows  from  the  equations  of  motion  of  an  ion  with  initial 
coordinates  p^  -  .  One  finds 


sin  <f>0|  S  G(pq)/p0  , 


(21) 


G  =  [1  +  In  pj/a*  -  v  (pj  -  l)/a*pf]1/2 


m 


where 


a* 


m  IB * 


\p  6 
P 


_  f 


0* 


/  ■/ '  2 

ml(v^ 

2Z.kT 
1  e 


(22) 


For  >>  1,  0*  $  1,  as  assumed  here,  we  have  a*  <<  1.  For  all 

a*  <  1  there  is  a  value  p  such  that  G(p  )/p  ■  1  and  then 
G(p0)/pq  >  1  for  pQ  <  p^.  pA(^),the  boundary  of  A*  given  by  Eq. 
(21),  has  the  form  indicated  in  Fig.  2. 2  As  the  ion  temperature  in¬ 
creases,  a*  goes  up  and  point  q  moves  down  reaching  p  ■  1  at  a*  ■  1; 


for  a*  >  1,  we  have  G(pQ)/p0  $  1  for  all  pQ  *  l.  Now,  d\l>/dp  -v  p 


-1 


for,  say,  p  <  Pm/3  [see  Eq.  (16)]  and  the  potential  field  is  then 


logarithmic.  The  mean  velocity  of  an  ion  with  p^  in  that  region,  in 
its  trip  to  the  probe, is  nearly  independent  oi  p^;  specifically 
Pq/Tq  *  pm*  vrt'ere  tq(Pq»  *0*  *s  the  t*me  to  th»*  probe. 
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Therefore,  as  long  as  Pq  <  p^, 

a  2  2  “22 

a*  *  P«T  *  Owe  (23) 

In  p 

and 

j*/1  ^  eaVTipy2  a,  U)1/2t 6/e  .  (24) 

P  P 

j*/Jw  grows  linearly  with  t  until  t  -  =  tq(p0  -  p^,  -  ir/2),  when 

a*  switches  from  a  quadratic  growth  In  p  to  a  nearly  linear  one;  for 
even  larger  t,  j*/ja>  will  actually  decline  because  Ions  from  outside 
the  3iji/3p  ^  p  '  region  will  begin  to  be  collected  and  both  field  and 
ion  mean  velocity,  p0/tq,  will  rapidly  decrease.  Neglecting 
logarithmic  variations  we  get,  from  (21), 

Pq  *  (<i»p6/ 8*) 1/2  (25) 

so  that  putting  t  ■  t  'v-  p  /p  in  (24)  we  obtain  the  peak  in  1*/1  , 

q  .  q  m  *  «• 

I*/J.  *  (*p/e*)1/2  6  ,  (26) 

which  can  be  far  greater  than  unity  and  depends  on  0*.  The  overshoot 
represented  by  (26)  may  be  seen  as  caused  by  the  sudden  set-up  of  the 
potential  field  which  traps  low  angular  momentum  ions  in  the  neighbor¬ 
hood  of  the  probe;  when  6*  decreases,  pq  increases  and  the  low  angular 
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momentum  region  Increases  too.  Eventually,  however,  a  value  of  0*  Is 

reached  for  which  moves  out  of  the  3i|i/3p  ^  p  *  region;  the  growth 

of  j*/jco»  38  given  in  (24),  is  then  stopped  by  the  rapidly  decreasing 

value  of  Pq/Tq,  and  not  by  the  fact  that  the  boundary  of  A*  has  been 

reached.  The  effect  is  the  same  for  all  smaller  6*  so  that  the 

current  is  now  insensitive  to  the  actual  value  of  g*.  The  critical 

value  of  0*  is  found  bv  putting  p  'v  p  ,  i.e.  x  <v  1, 

q  m  q 


0*  ^  e 


the  maximum  current  peak  possible  is 


j*/l  ^  *1/2  6/ g  . 

00  P 


(27) 


(28) 


According  to  the  H-S  similarity,  the  functions  given  in  Eqs. 

(20)  and  (6)  (for  0*0)  are  the  same.  Equation  (20)  depends  on 
x,  g,  0  and  i|»p.  In  the  flow  problem  the  value  corresponding  to  t  is 

A  A 

x^  =  the  parameters  i  and  M  in  Eq.  (6)  appear,  therefore, 

combined  in  a  single  one.  All  the  qualitative  results  derived  above 
for  J/J  need  not  be  repeated  and  we  shall  or.ly  add  two  new  points. 
(1)  For  the  present  problem,  Eqs.  (19)  and  (20)  become 


N  e  U  a(i/U) 

J  2irr  1 

P 


(29) 


i_. 


- To"  «(*/m) 

(syi/z  i 


(30) 


-12- 


Equation  (19)  makes  clear  the  meaning  of  a(l/U):  The  total  current  to 
the  probe  is 

J  -  2irrpl  J  -  NQe  U  a(i/U)  (31) 

so  that  a(t/U)  is  an  effective  probe  cross  section  in  the  plane  per¬ 
pendicular  to  the  flow.  (2)  As  already  indicated,  the  analysis  of 
Ref.  6  is  only  valid  for  t  <  t  ;  one  finds  that  t  5  u  t  a*  3.  Thus 
the  results  of  the  present  paper  will  be  valid  for,  roughly, 


t  <  3M  .  (32) 

We  note  that  the  peak  in  j/j  occurs  at  t  %  1. 

The  theory  given  above  explains  the  large  value  of  j/j#  at  0  ■  0. 

It  can  also  explain  its  sharp  decrease  when  the  probe  is  turned  by  a 

small  angle.  For  0  +  0  the  problem  changes  in  three  respects.  First, 

the  unperturbed  distribution  function  in  the  p  -  <p  plane,  £q(^)»  has 

now  a  drift  velocity  U  sin  0;  second,  the  Laplacian  in  Poisson's 

-2  2  2 

equation  includes  the  term  p  3  iji/3$  ;  finally,  angular  momentum  is 

not  conserved.  </e  shall  nov  assume  that  the  last  two  changes  have  no 

substantial  effect  on  the  current  to  the  probe;  this  point  is  discussed 

in  the  Appendix.  Then  the  results  for  j*(v*)/j  ,  v*  arbitrary,  are 

A  •  * 

not  changed;  the  drift  is  taken  into  account  by  using  the  new  form  of 

<+  1/2 

f 0 (vj* ) ,  which  now  has  two  characteristic  velocities,  (xT^/m^)  and 
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U  sin  6.  It  Is  clear  therefore  that  Eq.  (6)  may  be  written  ai. 


i.  .  i_  ( - 1 - 

j  j  VM  cos  e’ 

J  00 


M2  sin2  6  . 

;i*  —2 - ’  V 


(33) 


Now  the  characteristic  values  of  0*  and  a*  to  be  used  in  our  earlier 

discussion  on  j*/j  are  not  8*  ^  0  and  a*  'u  20/i p  6  but,  say, 

“  P 


2  2 
M  sin  6  _ 


6*  %  8  +  ~  ■  y- . .  i 


(34) 


2  .  M2  sin2  0. 

a*  %  — -  (8  + - 2 - ^  =  “t 

ip  6  ^ 

P 


(35) 


2  2 

(although  the  effects  of  &  and  M  sin  6/2  in  (33)  are  not  exactly 
additive,  they  may  be  considered  so  in  a  qualitative  discussion). 

As  long  as  6  is  so  small  that  By  %  8,  the  current  remains  fairly  con- 
stant.  When  0  becomes  of  the  order  of  (2B)  /M  or  (2e  )  '  /M, 
whatever  the  largest,  the  current  begins  to  decrease.  When  6  is  so 
large  that  aT  *  1,  we  have  Pq  S  1  and  so  G(pq)/Pq  in  Eq.  (21)  is 

A 

always  less  than  one.  Thus  a*  always  grows  (almost)  linearly  in  pQ, 
while  Pq/tq  remains  fairly  constant.  Therefore  %  1  (except  for 

£/M  cos  0  very  small) .  We  note  that  in  fact  the  peak  in  the  current 
disappears  for  6  smaller  than  the  value  for  which  aT  ■  1. 

Before  using  the  results  of  the  preceding  time-dependent  theory 
in  the  flow  problem,  we  must  examine  more  carefully  the  assumptions 
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behind  the  ll-S  sincerity.  First,  if  the  82*/3c2  tern  is  retained  in 
Poisson's  equation,  Eq.  (13)  would  read 


HI 

3p 


-  »p  4(g) 


+  <v(p,o  -  r2-^-|>  (p2  _  i) 


3C 


(36) 


where  <  >  has  the  same  meaning  of  Eq.  (14)  and  we  have  used  the 

variable  C  equivalent  to  x.  We  want  to  find  out  some  quantitative 
condition  for  writing 


<v  -  r2  32lj//aC2>  %  <v> 


(37) 


in  Eq.  (36).  The  equation  itself  can  be  used  to  obtain  32*/3C2  by 
integrating  once  with  respect  to  p  and  deriving  twice  with  respect 
to  5.  If  (37)  i8  valid  we  get 


a2^  6  ,2,-1  -1 

‘  3C2  ’  £  I2S  ln  0  -  2*<V>2> 


+ !  I  *  d°'  <&» 

1  Pm 


so  that 


<v-r24%y.i!(?(21no 

3C  2l2 


-  1)  < 


d26~1 

dC2 


26(^)2> 


+  °—i 
2om  3< 
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Both  (2  In  p  -  1)  and  p  have  maxima  at  the  largest  value  of  p  con¬ 
sidered;  thus  we  have 


**-i-*4~  **-*4 


3t2 


2dt2 


We  estimate  dd-1/dc  %  dd-2/ic  <  d~2  (£■»)-  i-1  (c  "  0), 

2-i  2  -1  2  2 

d  6  /dc  %  2A6  ,  3  v/3;  %  2Av.  For  the  typical  values  of  c  and 
i|»p  to  be  later  considered,  6  ^(5  ■  »)  -  6_1(?  *  0)  ^  0.15  and 
Av  <  0.20.  Thus  we  obtain 


3C2  PWm  t2 


Since  46  Is  never  far  from  unity  and  the  absolute  error  In  writing 

<v>  ■  v  is  tyoically  0.1  we  find  that  the  condition  for  the  neglect 
2  2 

of  3  iji/3c  in  Poisson's  equation  is  roughly 


l  >  * 


Second,  If  the  z-equatlon  of  motion  of  an  ion  is  twice  Integrated 
up  to  z  *  £.,  with  Initial  conditions  z  ®  0  and  dz/dt  *  U,  we  get  an 
equation  for  the  time  t^  that  the  ion  takes  to  travel  the  probe 


Ut*  Ut*  \ 

- -  f  1  +  - -  -a-1 

*  *  2M2 


where  i|t^  Is  given  by 
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(tt  -  t)  dt  3*/5>C 


(41) 


and  this  Integral  Is  along  the  trajectory  of  the  ion.  If  3i|>/3i;  Is 

small  enough,  the  bracket  In  Eq.  (40)  becomes  unity  and  we  get  the 

uniform  cotlon  assumption,  tfc  ■  i/U,  on  which  the  H-S  similarity 

_  o 

is  based.  This  can  only  be  true  if  is  very  small.  Integrat¬ 

ing  Eq.  (36)  with  respect  to  p  and  deriving  with  respect  to  c  we  get 


ip  In  p 
P 


.2  d6 


-1 


dC 


+  t 


(42) 


so  that,  at  most. 


%  0.4 
2M 


(43) 


where  we  considered  the  worst  possible  case  (initial  p  equal  to  p  ) 

m 

and  estimated  dd  */d£  Jfc  0.15,  3v/3c  %  0.20.  rrom  Eq.  (43)  we  con¬ 
clude  that  if 


<  M2  ,  (44) 

we  have  t^  ■  l/U  with  an  error  of  3 X  or  less. 

Third,  we  note  that  (1)  to  actually  obtain  the  current  to  the 

2 

probe  we  muct  add  to  Eq.  (31)  the  expression  N^e  Uirrp  [the  front  end 
of  the  probe  was  excluded  from  a  (l/U)],  and  (2)  Sanmartin's  analysis 


ffetJuniKM 
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l,: 

concerned  a  probe  whose  potential  was  switched  from  zero  to  V  at 

P 

t  ■  0  (Instead  of  a  probe  suddenly  immersed  in  a  plasma).  The 
corrections  resulting  from  points  (1)  and  (2)  are  only  important  if 
t.  <<  1,  but  will  be  incorporated  into  the  results  in  the  next 
section. 

Ill  -  The  Ion  Current 

An  expression  for  j^/j^  was  obtained  in  Ref.  6.  Neglecting 
some  terms  that  amount  to  less  than  2%  and  taking  into  account  the 
corrections  indicated  in  the  last  paragraph  of  Sec.  II,  we  get 
(for  t  -  t£) 


J- 


<8V1/2 


ir  a 


a*  <  h 


(45) 


(®V  Tt 


+  —  9  — - - 

a*  ,  : 

In  o 


_~2 


2  ,  -l,h  .1/2 

{1  -  -  sin  (-j-) 
it  a" 


f(ir“  d1/2}] 


a*  >  h  ,  (46) 


where  a  is  given  by 

a  Erf (In  a)1/2  -  (2*  7)1/2  F(t5)/c 

P  *■ 


(47) 


and 


-18- 


h  -  In  o2/(o2  -  1) 


(48) 


*<Tt)  -  Tjt'ir  +  0.6  t2)~1/2 


(49) 


The  actual  nondimensf onal  Ion  currert  j/j  may  be  obtained  from 


j/l 


I  fo(^! 


]•/]. 


(50) 


The  unperturbed  ion  distribution  function  in  the  plane  z  *  0  is 
Mav.ellian  with  a  drift  U  sin  0  so  that  from  Eq.  (50)  we  get 


2  71  oo 


t'l  dY  I  V*  dv?  2^ 

n  rv  A 


-  ra. 


0  0 


exp  1YktJ  {(1*)2  +  1,2  sin2  0 


-  U  sin  0  cos  y)J  j*/J, 


V  'xp  <-  ?  “p  <-  £->  'o'2  £  vl/2)  f 


(51) 


where 


n  *  2B/i i>  5  h  ,  y2  -  M2  sin2  0/ij»  5  h 

P  P 


(52) 


Using  Eqs.  (45)  and  (46)  we  can  write  (51)  as 
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*  TTTT1^ —  Ix°(u'  n)  +  3Xl<1‘'  n)1 

"  <8*p)  Tt 


(53) 


where 


s  ■  2(1  -  a  2)/ln  a 2 


and 


xQ(y,  n)  =  -  n  1  exp(-  y2/n) 


1/2 

dv  exp(-  v/n)  IQ(2yv  '  /n) 


+  n 


-1  <*p(-  -L)  | 


1/2  .  .  .  ,  ..1/2 
dv  exp(-  X)  IQ(2y  *j->  -  [sin-1  — /T  +  ^ - ]  , 


(54) 


X^y,  n)  =  xjn  -  n  1  exp(-  ^-)  j 


dv  exp(-  I  (2yv1^2/n) 

n  u 


1  ..  2,-1  1  2  ,  ,.1/2, 

•  7  11  -  7  sln  172  ■  7  <v  -  u  1 

V 


(55) 


It  has  not  been  possible  to  carry  out  the  Integrations  In  (54) 
and  (55)  analytical ’y.  A  number  of  limiting  expressions  may  be 
easily  derived,  however.  For  n  ■*>  0  (cold  ions  limit)  we  get 


xJJ°  +  sx“°  -  1  , 


V  (  1 


i  ,.„-i  i  +  hLls£!1}  + 1_ 

IT  y  y  2 


(56) 


n  2  ,  -1  1  2,2  , N 1/2, 

11  “7sin  v  ‘  »  (w  '  1}  1 


y  >,  1  • 
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For  11  ■*  0, 


X®*  +  sxjn 


Erf  n“1/2  +  2(ffn)"1/2  exp(-  n"1) 


2n“1  Erfc  n"1/2 


+  sin"1  E  1(2n“1)  -  2(irn)"1/2  axp(-  n"1)  Erfc  n'1/2]  .(57) 


For  n  fixed  and  u  ■+•  • 


Ypn  4.  a.  Y^0  X  »vW°  «  4  2 

Xn  +  aX.  %  X.  +  aX.  % - s  —  ; 

0  1  0  1  ^  iry  wp  * 


(58) 


all  curves  approach  the  cold  ion  limit.  Finally  for  p2/n  fixed  and 


n  * 


yHH  .  aYun  a.  4 

x0  sXl  %  ,1/2 

(im) 


«P<Z2^»  Vfr  (1  -f>  • 


(59) 


2  —1/2 
as  p  /2n  -*•  0  this  equation  approaches  4(irn)  (1  -  s/2),  which  is 

the  limit  of  Eq.  (57)  as  n  "•  Equation  (59)  can  be  also  rewritten 

as 

C  +  «r  %  if  4,1/2  <2')1/2  a  -  f)  :  (60) 


as  p  /2n  Eq.  (60)  approaches  4(1  -  s/2)/wp,  which  is  the  limit 


of  Eq.  (56)  as  p  ■»  ». 
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X  =  Xq  +  sX^  Is  given  graphically  in  Fig.  3  as  a  function  of  p 
for  several  values  of  n.  For  each  n,  curves  for  two  values  of  s 
have  been  represented  (s  ■  0,  s  ■  1/3);  interpolation  and  extra¬ 
polation  for  different  s  are  Immediate  because  X  is  linear  in  s. 

We  note  that  X  is  practically  always  very  close  to  X^.  The  function 
X(p)  is  a  direct  representation  of  the  peak  structure,  since 
X  'v*  j/j^  ind  p  ^  0. 

In  the  computation  of  X  an  overshoot  was  observed  for  the  largest 
values  of  n  in  Fig.  3:  in  approaching  the  cold  ions  (n  ■  0)  curve, 
each  (large)  n  *  constant  curve  overshot  it  and  then  approached  it 
from  above.  This  effect  was  so  small  that  it  could  not  ^how  up 
clearly  in  the  figure,  and  all  curves  were  interrupted  when  first  meet¬ 
ing  the  n  *  0  curve.  The  existence  of  the  overshoot  may  be  seen 
explicitly  in  Eq.  (60),  valid  for  large  n,  since  the  function 

(2ny)1/2  exp(-  y)  IQ(y)  , 

which  is  zero  at  y  ■  0  and  unity  at  y  ■  »,  has  a  maximum  1.17,  at 
y  ^  0.80.  Actually,  this  would  indicate  that  the  overshoot  should  be 
substantial;  thus,  one  may  conclude  that  for  the  moderately  large 
values  of  n  here  considered,  finite  n  effects  partially  mask  the  over¬ 
shoot. 

The  fact  that  the  n  ■  0  curve  is  not  an  upper  bound  of  the  family 


n  ■  constant  implies  that  the  seemingly  obvious  condition  3X/9ri  <  0 
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is  violated  for  some  values  of  n  and  p.  lh.it  this  is  possible  may  be 

easily  understood  by  noticing  that  if  v*  is  the  vectorial  composition 

of  an  equiprobably  oriented  (thermal)  velocity  and  a  directed  velocity 

(the  drift)  forming  an  angle  y  with  each  other,  an  increment  of 

either  component  may  decrease  the  value  of  v*  for  a  certain  range 

of  values  of  y;  under  some  conditions  the  increase  in  current  in 

this  range  may  dominate  the  decrease  that  appears  at  all  other  values 
of  y. 

From  Eqs.  (56)  and  (57),  closed  formulae  may  be  derived  for  the 

main  features  of  the  current  peak,  that  is,  its  maximum  and  its 

* 

angular  half-width  6^^  (the  width  of  the  peak  at  half-value  of  its 
maximum) .  For  the  maximum  we  have 


i_ 


max 


4~  (0  *  0) 

J  m 


ireo 


<%> 


1/2 


[X, 


On 


BXJ"] 


(61) 


the  functions  X^11  and  X^n,  given  in  Eq.  (57),  are  graphically 
represented  in  Fig.  4  for  convenience.  For  the  half-width  we  have 
the  condition 


+  sX 


un 


0.5[X°n  +  sxjn] 


where,  from  Eq.  (52), 


P 


M8i/2(o2  -  1)1/2 

2(*  1  In  a2)112 

P 


(62) 


(63) 
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We  now  note  that  all  curves  in  Fig.  3  meet  the  cold-ions  limiting 

a 

curve  at  values  of  p  clearly  smaller  than  p^^*  Thus  we  can  rewrite 
Eq.  (62)  as 


Xq0  +  sXj°  -  0.5  [X°n  +  sxjn]  ; 


(64) 


this  equation  only  involves  the  functions  given  in  Eqs.  (56)  and 
(57).  A  useful,  explicit  approximation  for  maX  be  obtained  by 

neglecting  the  dependence  on  s  and  writing 


xjj°  %  4/ tt  [p2  +  (4/tt)2  -  1]"1/2 


which  has  an  error  of  less  than  3%;  we  then  have 


ip  6  In  o  .  _ 

e, %  2I-V-5 - ]  [(  8 


'1/2 


.M2(o2  -  1) 


"X, 


ri2  +  1  -  £>211/2 


(65) 


Equation  (64)  has  been  solved  exactly  for  s  *  0  and  s  *  1/3; 
is  given  graphically  in  Fig.  5. 


Figure  6  presents  j/j 


versus  i/M  for  a  fixed  ip  and  several 
max  Tp 


values  of  e,  from  both  Eq.  (61)  and  the  experimental  data  discussed 

_2 

in  Ref.  6.  Theoretical  curves  are  presented  for  both  6  -  10  (full 

_3 

line)  and  8  a  10  (dashed  line).  In  Ref.  6  it  was  estimated  that 

_2 

in  the  experiments  6  was  of  order  10  or  less,  and  it  was  assumed 
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Chat  a  cold  Ion  theory  would,  therefore,  apply.  Our  analysis  shows 

that  the  condition  for  a  cold  ion  theory  is  not  B  <<  1  but  n  <<  1  (or 

more  weakly,  n  <  0.5,  say);  this  shows  up  clearly  in  Fig.  6  for  the 

largest  values  of  e  ^  and  l/M.  The  agreement  with  the  experiments 

_2 

is  excellent  for  e  -  0.009  and  0.041  if  6  10  ;  for  e  -  0.08,  the 

error  is  no  more  than  20%  (except  for  a  datum  obviously  in  error), 

_2 

still  within  the  error  of  the  measurements.  If  8  were  2  x  10  , 

say,  the  overall  agreement  would  greatly  Improve. 

Hester  and  Sonin's  experiments  exhibited  a  linear  dependence  of 

the  current  on  the  potential;  this  is  also  in  agreement  with  our 

1/2  —  1/2 

theory  since  jm  *»*  and  j/j^  .On  the  other  hand,  Bettinger 

—  3/2 

and  Chen's  theory  predicted  j  ^ 

Figure  7  presents  a  nondimens ional  half-width  versus  fc/M  for 

-2 

the  same  conditions  of  Fig.  6,  from  both  theory  (fi  line,  8  ■  10  , 

-3 

dashed  line,  8  ■  10  )  rod  experiments.  The  half-width  0^^  *8  not 

A  A 

ei/2;  it  is  defined  in  the  same  way  of  0jy2*  except  that  now  the  peak 
is  defined  as  the  current  in  excess  of  that  predicted  for  an  infinite 
probe.  1®2./2^BC  is  the  Pre^iction  Bettinger  and  Chen's  theory, 

for  cold  ions.  The  use  of  8j./2  allows  direct  comparison  with 
the  experiments  of  Ref.  6.  The  agreement  is,  in  general,  good  for 
8  *  10"2. 

The  experiments  showed  no  dependence  of  ®]./2  ^or  °l/2^  °n  ^p’ 
this  is  also  in  complete  agreement  with  our  theory,  as  easily 
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verified  in  Eq.  (65)  (on  the  other  hand  the  theory  of  Ref.  2  predicted 
-1/4 

0l/2  %  a^so  note  that  for  cold  ions,  our  theory  predicts 

that  8^/2  *-8  linear  in  e;  this  would  explain  the  claim  in  Ref.  6  that 
the  experimental  data  for  ®jy2^0l/2^BC  correiate  in  a  single  universal 
curve.  However  we  point  out  that  for  many  of  the  experimental  points 
of  Fig.  7,  as  in  Fig.  6,  the  ions  can  not  be  considered  as  cold 
(n  >  0.5).  [Even  for  cold  ions,  all  curves  do  not  exactly  meet 

A 

except  for  small  i/M,  because  of  logarithmic  effects]. 

IV  -  Conclusions 

The  present  paper  deals  with  a  significant  end  effect  in  the 
current  response  of  a  cylindrical  Langmuir  probe  in  a  collisionless 
plasma  flow.  Inf initely-long  probe  theory  predicts  that  when  the 
angle  0  between  probe  axis  and  flow  direction  decreases,  the  current 
experiences  a  smooth  decrease;  for  a  finite  probe,  however,  the 
current  may  exhibit  a  strikingly  different  behavior,  in  the  form  of  a 
strong  peak  at  small  0.  The  peak,  which  may  be  substantial  even  for 
very  long  probes,  appears  when  the  potential  is  highly  negative  and 
both  the  ion-acoustic  Mach  number,  M,  and  the  ratio  of  Debye  length 


to  probe  radius  e  are  large. 


The  only  analysis  of  this  end  effect  available  until  now  were  a 

A  A 

crude  theory  for  the  regime  i  >  3M  (£  being  the  ratio  of  probe  length 
3 

to  Debye  length)  and  some  numerical  computations  for  £.  <  3M  and 
0  -  0.5  Here,  the  regime  E  <  3M  is  rigorously  studied  for  0  arbitrary 


} 
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wit  h  in  the  peak  region.  The  unperturbed  ion  distribution  function  is 
supposed  to  be  Maxwellian.  It  is  found  that  the  nondimens ional, 
average  current  density  j/j  [where  J  -  J(i  -*•  »))  may  be  written  as 


■  4~  Irr.  e,  6,  *  ,  M9] 
J.  M  P 


(66) 


where  8  is  the  temperature  ratio  and  -  -  eV  /kT  (V  being  the 

P  pep 

probe  potential  and  the  electron  temperature).  Specifically  it 


is  found  that 


(8 O172  * 


2  ]  *  X(y,  n,  s) 

“  J 


itcM  a 


(67) 


where 


v  -  [ 


(O2  -  1)1/2  M 
(*p  J  In  o2)172 


1  e 


(68) 


,  2J.Cq.-H 

n  -  2  * 
i|»p  6  In  o 


(69) 


s  -  2 


l-o 
In  o^ 


-2 


(70) 


and 


6  1  ■  In  e  1  +  Y(t,  y  , 


(71) 


a  Erf  (In  o) 


1/2  <2*P  4,1/2  */M 


c[ir  +  0.6  (£/M)2}172 


(72) 
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Y  is  given  in  Fig.  1.  X  is  linear  in  s  so  that  a  single  graph  for 
X(u)  with  two  families  of  curves  covers  all  conditions;  this  graph 
is  given  in  Fig.  3.  The  brackets  in  Eqs.  (67)  and  (68)  are  scaling 
factors  so  that  Fig.  3  gives  £J/jJ(0)  directly.  Explicit  formulae 
are  presented  for  the  main  features  of  the  peak  (its  maximum  and  its 
half-width).  All  predictions  of  the  theory  agree  well  with  experi¬ 
mental  data. 

The  end  effect  may  be  advantageously  used  for  diagnostic  purposes. 

The  first  point  to  note  is  that  the  electron  temperature  has  no  effect 

on  the  peak;  the  dependence  of  j/j  on  T  cancels  out  except  by  way 

00  0 

of  6;  such  dependence  is  very  weak  since  the  logarithmic  variation  in 
the  term  In  e  ^  is  substantially  balanced  by  Y,  Since  does  not 
depend  on  Tg  either,  other  aspects  oi  probe  response  (usually  the  slope 
of  the  logarithm  of  the  current  for  weakly  negative  potentials)  must 
be  used  to  determine  the  electron  temperature. 

Apart  from  parameters  not  related  to  the  plasma  (!i,  i,  rp)  and 
weak  (logarithmic)  effects,  the  maximum  and  the  half-width  of  the  peak 
can  be  written  as 


L. 


max 


1_ 


max 


3 


ZJ  eV 

_i _ E 

m . 


m. 


(73) 


(74) 


where  the  right-hand  sides  are  explicitly  known  functions* 


Nq  is  the 
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plasma  density  and  Z^,  m^  and  T^  are  the  charge  number,  mass  and  tem¬ 
perature  of  the  ions.  If  ja>(0  ■  0)  can  be  accurately  extrapolated 
f'-om  experimental  data  at  moderately  small  6,  then  j^Cn/2)  and 
CO)  jointly  yield  TJ  and  Z  eV  /m  .  Equations  (73)  and  (74)  then 

00  0  l  p  X 

give 


VP-  W  Vzi 


2 

If  that  extrapolation  is  not  possible  or  2<|>  /M  is  large  [j  (0)  %•  j  (ir/2)], 

p  00  u  00 

some  additional  datum  is  necessary  to  find  all  four  N^,  V  ,  rn^/Z^ 

and  Tj/Z^.  If  n  <  0.5,  however,  X  %  X(n  ■  0)  and  it  is  possible  to 

find  Nn,  V  and  m./Z.  without  determining  T./Z. . 

Up  i  i  i  i 

The  present  analysis  is  subject  to  a  number  of  restrictions. 

It  is  mainly  based  on  a  theory  recently  developed  for  a  closely 
related  problem  in  transient  probe  response.^  The  most  substantial 
restriction  on  the  validity  of  that  theory  is  that  one  must  have 


l  <  3M 


additional  restrictions  are 

e  «  1  ,  i|»  >>  1  ,  20/d;  7  <<  1 

P  P 

The  last  inequalities,  although  used  extensively  to  simplify  the 


mmnHBs* 
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analysis,  are  hardly  restrictive  however  because  the  end  effect  Is 
small  If  they  are  not  satisfied.  The  use  of  the  aforementioned 
transient  probe  theory  In  the  present  context  Is  found  to  also  re¬ 
quire  (approximately) 

*p  <  M2  - 

-2 

*p  <  1 

(Some  additional  conditions  are  discussed  In  the  Appendix.) 
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Append  lx 

When  the  probe  Is  not  exactly  aligned  with  the  flow  the  ion  density 
Is  not  centrally  symmetric  because  of  the  drift  U  sin  8.  In  Sec.  II 
It  was  assumed  that  this  asymmetry  had  no  sensible  effect  on  the  Ion 
current  to  the  probe.  The  asymmetry  shows  up  twice  In  the  equation 
for  the  radial  motion  of  an  Ion 


2  d^p 

:  dr2 


ii  +  £_  [A  +  E-2 

3p  3  lA0  +  e 


(Al) 


where  is  the  nond linens  Iona  1  angular  momentum  at  time  t  =  u^t  -  0 
(when  the  Ion  crosses  the  z  *  0  plane)  and  the  second  term  In  the 
bracket  repre  tents  the  change  In  angular  momentum  due  to  the  azimuthal 
field.  When  6  -  0,  we  have  *  0  and  it  is  possible  to  determine 

the  region  A*  and  to  find  the  function  a*(T)-in  Eq.  (20)  by  using 
Eq.  (16)  in  (Al) .  If  0  ^  0,  however,  we  have  3\J>/3$  +  0;  3iJ//3p  now 
depends  on  4>  and  the  second  term  in  the  bracket  In  (Al)  does  not 
vanish. 

To  get  an  estimate  of  the  Importance  of  the  asymmetry,  we  first 
note  that  in  the  plane  z  ■  0  the  ion  density  is  uniform.  As  z  in¬ 
creases  the  Ions  readjust  their  distribution  function.  If  the  probe 
is  sufficiently  long  the  "infinite"  (two-dimensional)  limit  density 
is  finally  reached.  For  simplicity  of  discussion  we  assume  that  the 
limit  charge  density  can  be  written  as 


v  -  v  +  v.  cos  & 

oo  c°®  9  00 


(A2) 
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where  v  and  v,  are  functions  of  p  only.  A  (very)  conservative 

c00  moo 

estimate  of  v,  for  p  <  p  would  be  v,  “1  [note  that  within  the 

p«°  in  poo 

0-range  Of  interest,  the  peak,  we  have  aT  <  1  or 
2 

(M  sin  0)  /2i b  <  6/2  <<  1];  for  p  >  p  we  may  take  v,  *  0.  If  the 

p  m  p00 

probe  is  sufficiently  short  the  asymmetry  has  no  time  to  develop. 

A  conserv;  tive  estimate  of  the  time  required  to  reach  the  limit 
charge  density  (A2)  would  be  an  ion  plasma  period  or  £  ■  2ttM.  The 

A 

longest  probes  here  considered  have  £  %  3M.  Thus  we  consider 
£  %  3M  and  assume 

v.(z  ■  0)  ■  0 
<P 

Vz “ l)  m  v cos  * h  *  v"1  * 

and  a  linear  growth  from  v,(z  ■  0)  to  v,(z  ■  £) .  It  is  then  possible 

p  p 

to  determine  the  <j>-dependent  part  of  i|>  and  its  relative  importance  in 
Eq.  (Al) .  We  find  that,  under  the  worst  conditions,  the  <|> -dependence 
affects  no  result  by  more  than  10%.  Because  of  the  conservative  con¬ 
ditions  used,  thxs  should  justify  our  neglecting  Bip/B<p  in  Sec.  II. 

There  is  an  assumption  underlying  the  preceding  discussion  and 
the  main  body  of  our  analysis  that  deserves  consideration.  This 
assumption  is  that  the  ion  distribution  function  at  z  *  0  is  the 
unperturbed  distribution  function  far  ahead  of  the  probe.  Obviously 
the  field  ahead  of  the  probe  tip  will  somehow  affect  the  ions  reaching 
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the  plane  z  *  0.  Since  e  <<  1  and  the  potential  field  around  the  tip 
should  be  roughly  spherically  symmetric,  we  ca.i  assume  this  field  to 
be  i>  %  ii>p/p,  (in  order  to  get  a  rough  estimate  of  the  importance  of 
that  effect).  Then  the  ion  motion  for  z  <  0  can  be  solved  exactly; 
we  find,  for  Instance,  that 


i  i 

1  (z  -  0)  -  1  + 

N0 

2W  *  (1  +  2<J>/M2)1/2 


1  +  2W 

-  1 


» 


It  may  be  shown  that  the  perturbations  in  ion  density  and  ion  azimuthal 

velocity  are  of  no  importance.  The  perturbation  in  ion  radial 

velocity,  however,  may  in  some  cases  he  so  large  as  to  affect  the 

ion  current  to  the  probe.  This  occurs  when  the  probe  is  very  short 
.  2 

(£  «  3M)  and  <Jip/M  is  close  to  unity  [remember  condition  (44)]. 

Finally,  since  our  entire  analysis  is  restricted  to  the  regime 

£  <  3M,  we  would  like  to  comment  briefly  on  Bettinger  and  Chen's 

~  3 

theory  for  the  regime  £  >  3M.  To  well  understand  point  (D)  below, 

A 

we  emphasize  here  that  the  limitation  £  <  3M  of  our  analysis  originates 

from  our  lack  of  a  good  approximation  for  3i|»/3p,  for  p  >  p  ;  such 

m 

A  ^ 

values  of  p  come  into  play  when  £  >  3M.  When  3iJ>/3p  is  unknown  a*(x), 
and  thus  j/j^,  cannot  be  determined. 

In  the  light  of  our  theory  and  of  Fig.  2,  Bettinger  and  Chen's 
approach  may  be  resumed  as  follows:  (A)  they  divided  region  A*  in 
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two  subregions,  call  them  and  AQ,  lying  inside  and  outside  the  sheath, 
respectively.  For  the  sheath  radius  [equivalent  to  our  or  exp 
(6  1)]  they  used  an  expression  "a"  patched  up  from  numerical  result's 
for  probes  in  quiescent  plasmas.  As  pointed  out  in  Ref.  5,  "a"  has 
a  wrong  dependence  on  (j^;  this  leads  to  the  wrong  dependence  of  j 


and 


©1/2  on  'p  indicated  in  Sec.  III.  (B)  They  assumed  that  all  ions 


in  are  collected;  this  leads  to  the  requirement  9  >  %  3M.  For 

^Q,  they  noticed  that  the  appropriate  equivalent  time  (probe  length) 

a 

was  not  t .  =  i/M  but  t.  -  t  .  (They  did  not  incorporate,  however, 

Jl  l  m 

this  correction  to  their  formulae;  instead,  in  a  comparison  to  some 

experiments,  they  adjusted  appropriately  the  experimental  data.  This 

may  lead  to  confusion:  In  the  resume  of  Bettinger  and  Chen's  theory 

given  in  Ref.  5,  the  correction  was  not  considered;  it  may  be  easily 

verified  that  with  such  correction  the  overall  agreement  in  Figs.  5a, 

b,  c  of  Ref.  5  would  improve).  (C)  To  compute  Tm»  power  laws  for 

Bip/dp  were  used  instead  of  the  correct  expression  given  by  Eq.  16. 

The  variation  of  t  with  the  particular  power  law  is  not  greatly 
m 

significant  however.  (D)  The  most  crucial  point  of  their  analysis  is 
that  for  a*  inside  they  wrote 


a*(T)  *  <Tt  -  Tm) 

where  the  proportionality  constant  is  obtained  from  the  limit  -*■  ®. 
This  is  obviously  wrong  because  the  value  of  at  the  sheath  boundary 
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is  still  0(1)  so  that  for  some  distance  outside  the  sheath,  ion  velocities 

are  not  well  approximated  by  their  asymptotic  values,  unless  6  »  1. 

(This  point  is  related  to  the  well  known  Bohm's  sheath  criterion.)  For 

8  <<  1  the  error  can  be  substantial  except  if  t.  is  very  large.  (E)  A 

final  error  is  the  approximation  of  p^($)  everywhere  by  its  asymptotic 

value.  The  approximation  is  good  for  Aq  but  not  fcr  A^.  Thus,  the 

resulting  error  should  be  noticeable  for  moderate  t^. 

Nevertheless,  on  the  whole,  the  agreement  of  the  theory  of  Ref.  3 

with  experiment  is  good  for  moderate  potentials.  Thus,  it  should  be 

possible  to  use  that  theory  f6r  l  >  3M,  provided  that  an  appropriate 

shsath  radius  (p  )  be  used  and  the  correction  t.  ■+•  t„  -  T  be  in- 
m  it  m 

corporated  into  the  formulae. 
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Fig.  2  Diagram  showing  regions  A*  and  B*  for  simplest  unperturbed  ion  distribution.  A* 
is  the  set  of  the  positions  at  z  »  0  of  those  ions  that  will  strike  the  probe  at 
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Fig.  6  Comparison  of  maximum  current  density  (0  *  0)  from 

experiments  (Ref.  5)  and  present  theory  (full  lines 

-2  -3 

6  ■  10  ,  dashed  lines  6  ■  10  ). 


experiment: 


Fig.  7  Comparison  of  angular  half -width  (as  defined  in  Ref.  5)  from  experiments 

-2  -3 

(Ref.  5)  and  present  theory  (full  lines  3  ■  10  ,  dashed  lines  B  ■  10  ) 


